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Abstract 

We present a study of planar physical solutions to the Lorentz-Dirac equation in a constant electromagnetic 
field. In this case, we reduced the Lorentz-Dirac equation to the one second order differential equation. We 
obtained the asymptotics of physical solutions to this equation at large proper times. It turns out that, in 
the crossed constant uniform electromagnetic field with vanishing invariants, a charged particle goes to a 
universal regime at large times. We found the ratio of momentum components which tends to a constant 
determined only by the external field. This effect is essentially due to a radiation reaction. There is not 
such an effect for the Lorentz equation in this field. 
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I. INTRODUCTION 



The Lorentz-Dirac (LD) equation has an ill fame to suffer from the various type inconsistencies. 
The latter come from the higher derivative Schott term entering the LD equation and appear as 
the blowing up (runaway) and acausal solutions. However, despite of its undesirable features, we 
have to accept the LD equation as a correct one in its range of applicability by the following 
reasons. First, as it was shown in the seminal paper by Dirac and then more elaborately in [Q], 
the LD equation stems from the energy-momentum conservation law provided a charged particle 
is sufficiently small and possesses neglible higher multipoles of a charge distribution. Second, 
under these assumptions the LD equation is a minimal evolutionary equation describing a radiation 
reaction which complies with all the symmetries of the model: the Poincare and reparameterization 
invariance. Furthermore, having made certain approximations, the LD equation was derived in the 
context of quantum electrodynamics (see, e.g., Q), where the LD equation can be considered as 
a leading quasiclassical asymptotics to the Schwinger-Dyson equations for an electron. Therefore, 
the LD equation makes a physical sense and, under certain conditions, its solutions should give rise 
to predictions which can be observed in experiments. It is clear that the LD equation is valid in 
the range of energies and field strengths where the quantum corrections are neglible in comparison 
with the classical contribution. Rough general estimates of this range can be found, e.g., in [||. ||], 
and the more accurate analysis for the particular case of a constant homogeneous magnetic field is 
presented in Various generalizations of the LD equation to include spin and higher multipoles 
jjj], or an interaction with non-Abelian gauge fields || and gravity to higher dimensions Jlpfl , 
to dyonic [11 1 and massless charged particles [ [12[ | are also known. All of them have the higher 



derivative terms and, hence, possess the same unwanted properties as the LD equation. 

There is a coherent approach [[|, p~3| P~7|] how to extract a physical information from the LD 
equation and its analogues. It is based on the notion of a physical solution. In a general setting, it 

looks as follows. Given a system of interacting fields 4>Z with the action functional 

N 

s[f x , . . . , <jP N ] = K's^Z) + s int [4>l ■■■All (i) 

a=l 

where 7 is a condensed index representing group and spacetime indices and spacetime points, a 
numerates the fields, and A a are some constants. Then the solution </>2(A) of the coupled system 
of equations of motion corresponding to the action (||) at given initial and boundary conditions is 
called physical if there exist finite limits 



lim^(A), a = l,N, (2) 



other lambdas being fixed. This regularity condition completely rules out the runaway solutions 
to the LD equation. Besides, its physical solutions are unambiguously specified by the six initial 
data - three position coordinates and three momentum components - as it should be in the realm 
of Newtonian mechanics. 

Since the LD equation is nonlinear, it is hard to solve it even in simple external field configura- 
tions. Almost all the exact solutions to the LD equation can be found in [0| In this 



paper, we address the problem of a finding and description of exact physical solutions to the LD 
equation in a constant homogeneous electromagnetic field. Moreover, we restrict ourself by a planar 
motion only. Even in this rather simple situation we did not succeed in finding the exact essentially 
planar (i.e. non linear) solutions to the LD equation. However, we reduce the LD equation to the 
one second order differential equation and investigate the asymptotics of its physical solutions at 



large times. For the constant homogeneous magnetic field this asymptotics was found in [15|. As 
far as the constant electric and crossed fields are concerned these asymptotics, to our knowledge, 
are obtained for the first time. It turns out that, in the crossed field configuration, the LD equation 
possesses an attractor and the system passes into a universal regime at large times. After the lapse 
of time, the identical charged particles moving on the plane in such an electromagnetic field "forget" 
their initial data. Their trajectories become parallel and a certain ratio of momenta components 
tends to a constant which is independent of the initial conditions and determined by the external 
field only. This effect is essentially due to the radiation reaction. It is absent for the solution to the 
Lorentz equation in this field and can serve as an explicit manifestation of a validity of the classical 
radiation reaction theory in the domain of its applicability. 

The paper is organized as follows. In Sec. O, we present general formulas regarding a radiation 
reaction and define the physical solutions to the LD equation. Here we also give an integro- 



differential equation for the physical solutions. Sec. [II is a main part of the article. In Sec. Ill A 



we briefly describe the linear motion and exact solutions to the LD equation in this case. In Sec. 



Ill B 1, we investigate the symmetries of the LD equation and provide the necessary and sufficient 



condition for the motion of a charged particle to be planar. In Sec. Ill B 2, we derive the second 



order differential equation describing planar solutions to the LD equation. Sec. I1IB2 is devoted 
to the asymptotics of the physical solutions to the LD equation at large times. In Sec. I [IB 4, 
we consider the same problem in the framework of the so-called Landau-Lifshitz equation [20]. In 
conclusion, we summarize the main results of the paper and discuss the prospects for a further 
research. 
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II. GENERAL FORMULAS 



Consider a particle with the charge e and mass m interacting with the electromagnetic field 
on Minkowski background M 1 ' 3 with the metric r]^ u = diag(l, —1, —1, —1) and coordinates x^, 
fi = 0, 3. The action functional for such a system has the form 



S[x(t),A(x)} = -m J drVx^-e J drA^x" - ^- J d 4 xF /IU F fJ-u , 



(3) 



where x^(t) defines the particle worldline, F^ u := d^A^ is the strength tensor of the electromag- 
netic field (the square brackets denote an antisymmetrization without 1/2) 



F, 



/it/ 



E x E n 



E 7 



-E x 



-Ey H z 



-E z —H y H x 



H Z Hy 



-H x 




(4) 



and we take the system of units in which the speed of light c = 1. In the proper time parameteri- 
zation x 2 = 1, the LD equation Jj], ^l] reads as 

2 



(5) 



where is the strength tensor of the external electromagnetic field. Introducing the dimensionless 
quantities 



x^ ->■ m r — )• m 1 r, — >■ m 2 e 1 F At!/ , 



we rewrite it in the form 



J/i = fv + ^v^ + v v,j,), U-=F liV v v , 



(6) 



(7) 



where A := 2e 2 /3 and mv^ := mi' 1 is the 4-momentum of the particle. 



The LD equation possesses unphysical solutions. Following [13, 14], we shall call the solution 
x^(A, t) physical if it tends to the solution x^(0, r) of the corresponding Lorentz equation as A goes 
to zero. It is a realization of the general definition given in Introduction in the case of classical 
electrodynamics. According to this definition, the physical solution should be regular at the large 
mass m and small e 2 . From Eq. (^) we see that this requirement leads to a regularity of the 
physical solution with respect to the external field and the parameter A entering (0). The former 
simply follows from the general theorems regarding dependence of solutions to ordinary differential 
equations on a parameter (see, e.g., p2| ), while the latter condition is not trivial. Also notice that 



all the known physically reasonable solutions to the LD equation are physical in the sense adopted 
by us. Some extra arguments in favor of this definition of physical solutions are given in Appendix. 

We can find these solutions perturbatively as a (formal) series in A. This perturbative scheme 
reduces the order of the LD equation and provides a unique solution to it at some fixed initial 
position and velocity of the particle. The first iteration of this perturbative procedure yields the 
Landau-Lifshitz equation |20| |. It is not difficult to write an integro-differential equation which 
describes the physical solutions to the LD equation |[|, p~7f] . In the proper time parameterization, 
it looks like 

/•oo 

v^t) = pr^(r) / dte^Vuir + Xt), pr^ := ^ - (8) 
Jo 

where = + XiP'V^,. Solutions to Eq. (^) are solutions to the LD equation (0) with v^v^ = 0. 
It is the latter requirement that gives rise to the projector entering Eq. (||). The solutions of this 
equation tend to the solutions of the Lorentz equation at A — > 0. Expanding Eq. (0) in a series in 
A, we see that solutions to Eq. (||) are those solutions to the LD equation which are obtained from 
it by the aforementioned perturbative scheme. If we knew all the terms of the perturbation series 
for the acceleration v(t, A) then formula (^) would tell us that this series in A must be summed by 
the Borel method pq| . So, if the following conditions are satisfied at some fixed initial position and 
velocity 

1. There exists a unique solution to the corresponding Lorentz equation, which is defined at 
any r > To and tends to infinity not faster than Me aT at r — > oo. Here, To, M and a > are 
some constants; 

2. The perturbative series in A converges absolutely in a vicinity of the point A = at sufficiently 
small A; 

then a unique solution to Eq. (||) exists at sufficiently small A = Ao > 0. If the value of Ao is smaller 
than the physical value of A then the physical solution to the LD equation at the physical value of 
A is obtained by an analytical continuation in A. 

A concrete prescription how to construct this analytical continuation depends on analytical 
properties of the solution to the Lorentz equation. For example, if this solution satisfies the first 
condition above and has a finite number of singularities in the part of the complex T-plane where 
ReT > To and IniT > 0, or ReT > to and Imr < for some To, then x^(A,t) has the same 
properties at sufficiently small A. In that case, we can rotate the ray along which the integration 
contour tends to infinity so as to make the integral (||) convergent for any a. In particular, this 

5 



procedure makes the right hand side of Eq. (|8|) convergent when we perturbatively solve Eq. (||) 
by Picard iterations starting with the solution to the Lorentz equation x^(0,r) satisfying the first 
condition above, while 3aX > 1. Of course, the latter situation is rather unphysical since 

x H m? 
aA~a— , H = — , (9) 
Ho \e\n 

where H is a characteristic value of the field strength, a is the fine structure constant and Hq is 
the Schwinger field. However, we can define the physical solution even in this case. 
Another, possibly more convenient, form of Eq. (g) can be derived if we write 

roo 

v 2 = dte-\vf){T + \t/2), (10) 



for physical solutions. Then 



^m( t ) = P r u( T ) / dte 



DC 



(11) 



f v (r + kt) + X ds(vf)(T + Xs/2)v u (T + X(t-s)) . 
Jo 

In this form, it is obvious that, in the absence of external fields, physical solutions to the LD 
equation are straight lines in the spacetime. 

The integro-differential equations for physical solutions to the LD equation, which we have 
presented in this section, are not very useful in finding analytical solutions. But they are pertinent 
to numerical simulations, for example, by the Picard iterations, rather than the LD equation itself. 
Even if we set the initial conditions to their physical values and would solve numerically the LD 
equation then we shall obtain an unphysical runaway solution owing to machine approximation 
errors. Also, Eqs. (B) and ( |Tl~D are a good starting point for an investigation of the stochastic LD 
equation (sec, e.g., ||, ^4|). As in the case of numerical simulations, the unphysical solutions should 
be explicitly excluded to get rid of stochastically induced runaways. The study of this problem will 
be given elsewhere. 



III. PHYSICAL SOLUTIONS FOR A PLANAR MOTION 

In this section, we study the planar motion of a charged particle obeying the LD equation. We 
call the motion of a particle planar (linear) provided that the trajectory of this particle in the 
space can be made planar (linear) by an appropriate Lorentz transform. For the linear motion, 
the worldline of a particle lies in a two-dimensional plane of the spacetime. For the planar motion, 
it lies in a three-dimensional hyperplane. Throughout of this section, we mostly assume that the 
particle moves in the constant homogeneous external electromagnetic field, but some results can be 
generalized to a non-constant electromagnetic field of the special configuration. It will be mentioned 
in its place. 
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A. Linear motion 



Let us consider, at first, the hyperbolic motion, which is the motion with vanishing the LD force, 



Vfj, + v 2 v^ = =>• v 2 = —Lxj 2 = const 



v»(t) = cosh(wr) + ^ sinh(wT), a 2 = -/3 2 = 1, = 0. (12) 

The hyperbolic motion is the solution to the LD equation (|7|) with a constant and homogeneous 
external electromagnetic field if, and only if, 

&* = F^F p u v u = -wV. (13) 



By the use of canonical forms (see |25| and also below) of the tensor F^ u , it is not difficult to show 
that the last equality (the equation on eigenvectors) is fulfilled if, and only if, there exists a Lorentz 
frame in which the particle moves along 3- vectors of the electric and magnetic field strengths. In 
this system of coordinates, the problem reduces to a description of the linear motion. 
For a linear motion, the LD equation (|7|) is equivalent to (see, e.g., ||l], 15, [ll| ) 



, V = = X-j- , = = + E => v{t) = sinh(c 2 + Et + c x e x ' r ). (14) 
V 1 + v « T V 1 + v 

where v := x(r). A generalization of this solution to the case of the electric field depending on r 
is trivial. The solution (|Ti| ) with E = is also a general solution to the free LD equation written 



in the Lorentz frame, where the initial 3- velocity and 3-acceleration are parallel. The solution (14) 
becomes physical if we take c\ = 0. It is easy to verify that the solution (|l4]) satisfies the integro- 
differential equation (||) only at vanishing c\. In order to make the integral convergent at XE > 1, 
we have to rotate the integration contour in Eq. (0) as it was described in the previous section. 



B. Planar motion 

1. General consideration 



Let us turn to the planar motion. When one of the Poincare-invariants of the electromagnetic 
field is not zero, the strength tensor (|j) can be represented as 

= wiejf e^ 1 + W24 Me 3 > (e^) = rj a/3 , (15) 

where e^, a = 0, 3, is a tetrad of eigenvectors of the tensor (F 2 )u- The eigenvalue uj\ of the tensor 
(F 2 )y corresponds to the vectors e^, and the eigenvalue — uj\ corresponds to the vectors e^. In 
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terms of the Poincare-invariants of the electromagnetic field I\ = E 2 — H 2 and I2 = 2(EH), these 
omegas read as 



u\ = {yjq + Jf + h)/2, uo\ = (^/ 2 + If - h)/2. (16) 

The Lorentz transforms, which do not change the strength tensor (|l5|), constitute the group 
SO(l, 1) x SO(2). Its matrix representation is obvious. 

In the degenerate case, when I\ = I2 = 0, the strength tensor is given by 



F^=ueKl\ (e a e b ) 




-1 



(17) 



where e%, a = { — , 1}, are eigenvectors of the tensor {F 2 )y corresponding to zero eigenvalue. The 
normalized eigenvector is orthogonal to the vector eg, which is a normalized eigenvector of Ft, 
and (F 2 )„ corresponding to zero eigenvalue. These conditions determine the vectors e± and e 3 
uniquely up to adding the isotropic vector e_ and inversion. The factor oj can be included to e_, 
but we leave it in the expression so as to control the external field. The strength tensor ([l~7| ) is 
invariant with respect to the two-dimensional Abelian subgroup of the Lorentz group generated by 
the two elements 

Af = rf v + nefe} + |e!e^ Af = jf v + r 2 e [ ?e v } + ^ e - e -> ( 18 ) 

where r\ and T2 are the group parameters. Thus, in both degenerate and non-degenerate cases, we 
anticipate two integrals of motion of the LD equation provided the eigenvectors of the tensors F^ 
and (F 2 )v do not depend on a point of the spacetime. 

In order to obtain these integrals of motion, it is useful to introduce new variables adjusted to 
the action of the symmetry group. In the non-degenerate case, they are 

v =: y^coshV', vi =: ^/u^sinhip, \ 2 y/uh^'f, v 3 =: ^/u^cos 99, (19) 

where v a := (e a v), u e — Uh = 1, and (p and ip are the symmetry group parameters. Then, convolving 
the LD equation (0) with the eigenvectors, we have 

,00 df (20) 
u h {T)y{T) = C2e x ~ lT - J ^e-^-W^i). 

In the degenerate case, the analogous variables read as 

vi=:v_r 2 , v 3 =:v_ri, v+ =: vi 1 + v_(r 2 + r|), (21) 



where e+ is an isotropic vector orthogonal to and eg, and such that (e+e_) = 2. The respective 
integrals of motion become 

u(r)ri(r) = cie A ' T , 

/oo jj. (22) 
ye~ A C- r W(i), 



where n(r) := v_(t). The non- vanishing constants ci and C2 in Eqs. (20) and ( f22| ) correspond 
to unphysical solutions. So, they should be set to zero and we do not take them into account 
henceforth. 



Consider the particular case I2 = (see ||14| ). If I\ > then W2 is zero and we get 

v 2 = v 3 = or ip = 0. (23) 
The first case is a linear motion, the second case is a planar one. When I\ < 0, we have uj\ = and 

V> = 0. (24) 

In the degenerate case I\ = I2 = 0, 

ri = 0. (25) 

By definition of the variables ip and r\ as the group parameters, the particle can be confined to a 
plane by an appropriate Lorentz transform in these cases. 

Thus we have proved the following statement. If I2 = and the field strength tensor admits 
the representation (O) or ([H]) with the constant eigenvectors then a charged particle obeying 
the LD equation executes a planar motion. In a constant homogeneous electromagnetic field, the 
converse statement is also true. Namely, if a charged particle obeying the LD equation executes an 
essentially planar (i.e. non linear) motion then I2 = 0. 

2. Second order equation 

Now we investigate the planar motion in detail. Let us characterize this motion by a tetrad 

e&lr,<% = V*fl, ^ = 0, (e 3 ) 2 = -l, (26) 

where the indices a and j3 run the values 0,1,2. They are risen and lowered by the metric 
Vap = diag(l, — 1, — 1). The worldline of the particle and the external electromagnetic field ad- 
mit a representation 

afi _ (27) 



where £012 = 1 an d £ 2 = {±1,0}. The LD equation ([?]) is rewritten as (for the Lorentz equation 
see, e.g., g|) 

v a = oje^v^C + + v 2 v a ), v 2 = 1. (28) 

We see that the LD equation of a charged particle confined to a plane possesses a symmetry. This is 
a residue of the symmetry discussed above after the reduction to the plane. The residue symmetry 
group is constituted by the Lorentz transforms leaving the vector £ Q intact. So, if £ 2 < this 
symmetry group is isomorphic to SO(l, 1), and if £ 2 > it is isomorphic to SO(2). This symmetry 
allows us to reduce the problem of integration of the system of equations ( |28"D to an integration of 
the autonomous system of three first order equations or the one second order equation. 
To this end, we introduce new more convenient variables 

m a ._ £ aP7vp VT (29) 
In these variables, the LD equation (|28|) turns into a system of the first order equations 

Xm a = m a + w(£ Q — pv a ), v a = —£aP'y m ^v 1 , m a v a = 0, v 2 = 1, (30) 
where p := £, a v a . Then we introduce the invariants of the symmetry group action 
a = ^ a m a , b = p~ 1 e a ^' 1 ^ a mRV 1 = — 

(31) 

-2/ 2 t 2 2\ 2 a 2 

s = p (m — £ a ), u = p — £ . 

These invariants are dependent. From their definition, it is not difficult to obtain the identity 

a 2i±i!^ + 6 2 = _ sn _ (32) 

The invariants (|3l|) evolve according to the equations 



Ad = a — cou, u = — 2b(£ +u), 

1 4- f 2 n ■ \F 2 \ - 1 a 2 4- £ 2 h 2 

Xs = 2s(l + Xb) + 2wai-i^, = 6 + A^ 2 s + Aa 2 ^ = b - X . 

£ Z + U u u 



(33) 



The first equation in this system is the one of the equations of the system (2C) or (22) written in 
a differential form. The physical solutions are described by the one integro-differential equation on 
u(t). It is obtained from the above equations if we write 

a(r) = ^Wi), p\r)s{r) = T ^-'(H^i), (34) 



and substitute these expressions to the identity ( |32[ ) with the function b(r) taken from the second 
equation of the system (|33|). The solution of this integro-differential equation is specified by the 
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only one arbitrary constant. If this solution or some unphysical solution to the system (33) are 
known, we can integrate the LD equation. 

Indeed, in a general position the vector v a can be expressed as a linear combination of the 
vectors £ a , v a and £ a ^^v^ ■ The coefficients of this decomposition are certain functions of the 
invariants which are already known. Therefore, we need to integrate a system of linear equations 
with variable coefficients. So, 

V a = -eaR-^V 1 + —(£a -Pv a ). (35) 

u u 

Because of the orthogonality condition, only two equations are independent. Now we make a 



substitution to Eq. ( p5[) of the form (|l^), (21). In the case £ Q = (1,0,0), we have 

v o = p, v\ = y/u cos tp, V2 = \/usinip, (36) 

and 

uip = —a, (37) 
whence the momenta v a are found in quadratures. If = (0, 0, 1) then we substitute 

vq = \/ucoship, v\ = y/usinhip, v% = p. (38) 

This results in 

uip = a. (39) 

In the third case £ a = (1, 0, 1), we do 

v\ = rp, V- = p, pv + = 1 + ur 2 , (40) 

and arrive at 

ur = a. (41) 

The case of an arbitrary vector £ Q is reduced to the considered ones by a proper Lorentz transform 
of the tetrad indices. Notice that the equations of motion (|33|), (|37|), (|39|) and ( pi] ) are valid for a 
non-constant external field parameter oj. In accordance with our general considerations, physical 
solutions are specified by four constants - two constants specify the initial position on the plane, 
one determines u(j) and another one is needed to pick out the unique solution from Eqs. (|3~7|), (39) 
or ©. 
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Thus, we have to find the evolution of invariants described by Eqs. (33). In case of a constant 
lj, the autonomous system ( ^3|) is equivalent to the one second order differential equation on the 
function a{u) or its inverse u(a) 



a" 



\2au + £ 2 {a + um)]a' 2A 2 a 2 (£ 2 + u)a' 3 
2u(a — wu)(( 2 + u) u(a — uju) 2 

.. _ [2au + £ 2 (a + uu)]u 2 2A 2 a 2 (£ 2 + u) 

^ — ~7 7 w~Zo " \ ~T~ 



(42) 



2u(a — um)(£ 2 + u) u(a — uju) 2 

Even in the simplest case £ 2 = 0, when these equations can be cast into the form 

,, aa' 2a 2 a' 3 .. au 2 2a 2 ,. nS 

a = -. r - t 72 , u = — + -, (43) 

u(a — u) (a — uy u(a — u) (a — uy 

we have not succeeded in finding a general solution. 



3. Asymptotics 



However, we can investigate the asymptotics of exact physical solutions to the LD equation at 
large times. It is easily done in the coordinates where the vector field of the system (|33|) has no 
singularities. Making a change of variables 



a = ua, 



b = ub, 



we come to 



AS = a[l + 2Xb(f + it) 



Xb = b- X[a' z -b 2 (f + 2u)}, 



(44) 



-2bu{f + u). (45) 



It is not difficult to find the stationary points of this system. We are interested only in physical 
solutions and physical stationary points. A physical stationary point as a particular case of a 
physical solution should be regular in A. Again we have three cases. 

The case £ 2 = 1 is a planar motion in a constant homogeneous magnetic field [[l5|, 27]. The 
system (Ea) has two stationary points, one of them being physical 



a = -, 
9 



9-1 
2A • 



u = Q, g := 2~ 1/2 (1 + V1 + 16A 2 w 2 



1/2 



(46) 



Linearizing the system (45) in a vicinity of this point, we obtain the asymptotics of the exact 
solution to the LD equation 

5a = -u(0)Ae x ~ 1{1 - 9)T + e A ~ V 



. . , . 2ujt , , . „. 2ujt 
(ci + u(Q)A) cos h (c 2 + u(0)B) sin 



5b = -u(0)Be x l ( 1 -9)^ + e x ^ T 



[c 2 + u(Q)B) cos 



9 
2ut 

9 



(ci + u(0)A) sin 



9 
2ujt 

9 



(47) 



5u = u(0)e x 1( - 1 - 9)l 



A 



5(55-4)' 



B :-- 



3(ff ~ I) 2 
2A(5 5 - 4) 
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The terms in the brackets describe runaway solutions. They are unphysical and have to be set to 
zero by a proper choice of the initial conditions. Then the physical solution to Eqs. ( |45|) is solely 
specified by the initial condition on u. The first term in the first line in Eqs. ( [47| ) describes a 
correction to the rotational speed of a charged particle due to the radiation reaction. Inasmuch as 
u is non-negative, this correction has an opposite sign with respect to the main contribution, i.e., 



the rotational speed increases with time and tends exponentially to its limiting value (46). The 



limiting value is, of course, lesser than the cyclotron frequency. In case at hand, u is related to the 



kinetic energy of the particle and so the expression for 5u in (47) describes its decreasing. 

The case £ 2 = — 1 corresponds to a planar motion in a constant homogeneous electric field. The 
system ( f45| ) possesses two stationary points. Only one of these points is physical 

a = oo, b = ^—, u = l, g := v / l + 4A 2 w 2 . (48) 

The linearized in a neighbourhood of this point LD equation has the solution 

5d = -um(0)(l - g-^e^^ 1 "^ + e x ~' T [c x + wu(0)(l - g' 1 )) , 



5b = -u(0)Be x ^-fr + ^± e x ^ T \ Cl + w «(0)(l - g' 1 )] + e A V 

5-1 

5u = u(0y- 1 ^, ^ := (5-l) 2 (25 + l) 
V ; 2\g(2g - 1) 



2ciAa; u(0)(ff-l) 
° 2 5-1 X(l-2g) 



(49) 



Unphysical solutions are the terms in the square brackets. Demanding their vanishing, we uniquely 
determine the integration constants c\ and C2 through u(0). The correction to the "frequency" tp 
increases with time and tends exponentially to the limiting value (48). The expression for 5u(0) 



in (^) describes an evolution of the square of the momentum component normal to the electric 
field. As expected, this component exponentially tends to zero and the solution passes into the 
hyperbolic motion (12). 



The case is more involved. To shorten formulas, we redefine the variables entering (p5|) 

a — > wa, b — > Xoj 2 b, u — > (\oj)~ 2 u, r — > Ar, (50) 

and shall restore the original notation where it will be necessary to make estimations. After this 
redefinition, a regularity in A, which distinguishes physical solutions, means a regularity of the 
solution in t . Then the system ( f45| ) has a single stationary point 

5=1, 6 = 1, u = 0. (51) 
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This point is degenerate and, therefore, the solutions to the linearized system improperly describe a 
behaviour of solutions to the LD equation in a vicinity of this point. To obtain a correct asymptotics, 
we integrate the last equation in (|45|) 



u = u(Q) 



1 + 2n(0)r + 2u(0) / dtSb(t) 
o 



-i 



(52) 



The integrand of the third term in the square brackets tends to zero. Consequently, the second 
term in the square brackets will dominate at large times 



r>A, r > (2Au^u(0)) , 



(53) 



and we can take 



(54) 



u « t-' l /2. 

Then the equations for the leading asymptotics read as 

5a = Sa + T^db + r- 1 , 6l = -25a + 6b + r -1 , (55) 

where we keep only the leading terms. The system ( |55[ ) possesses runaway solutions, which are 
nonregular in A (t -1 ) and, consequently, unphysical. They can be removed by an appropriate 
choice of the initial data. As for the physical solutions, their leading asymptotics takes the form 

. .,_ > 3/r 

a = 1 - t- 1 - T- 2 (31nr + 2«(1) - 6) + o(t" 2 ) 



6=1- 3T- 1 
1 



T- 2 (91nr + 6u(l) - 23) + o(r~ 



t I 8u 3 (l) I 

3 / e- 2 °/V 
r I 8w 3 (l) 



3/t 



+ o(r" 



(56) 



atr" 1 + ^(Slnr + 2.(1) - 1)] + o^) = J- 



3/t 



+ o(r~ 



Here we also add a next to leading correction to the asymptotics, which can be derived from the 
initial non-linear system ( |45| ) . The last equalities in (^fj) show how the power of the proper time 
entering the asymptotics tends to its limiting value. The initial value u(l) in these last formulas 
differs from u(l) appearing in the second equalities in (j5q). These initial values are related in an 
evident manner. 

Substituting the asymptotics (56) into Eqs. (|40|), ( |4l| ) and bearing in mind the replacement 
(|5l|, we find 

T\A/2r 

-A/2r . s (57) 



'■'1 



T \l/2 

2A 



(2. 2 (A ) g A/2r + ((A/T ) V 2 ) 



= Ac (£) 



r \3/2 / Aw 2 r 



8*; 2 (A) 



+ ((r/A) 1 / 2 ) =-.; 2 . 
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As we see, the system goes to the universal regime. For example, the quantity 



v\ vf V Xv (X) 



-3---3-^[^7T^) (58) 



ceases to depend on the initial data and tends to Xu. This occurs on the proper time scales 



T > A 



In- 



(59) 



Xuivf(X) 

This effect is essentially due to the radiation reaction. After the lapse of a certain time, the charged 
particles moving on the plane in the electromagnetic field with the invariants 1% = I2 = will 
have the same ratio of momenta of the form (|58| ) irrespective of their initial momenta. So, if we 
measure the momenta components of these charged particles, the measured data will lie on the 
cubic parabola determined by Eq. (|58|). 

For comparison, we give here the well-known solution to the Lorentz equation in this electro- 
magnetic field 



vi = vi(0) + y/u(0)u)T, v = v (0) + wi(0)wt + >/u(0) 
v 2 = v 2 (0) - vi(0)ut - a/u(0) 



2 (60) 

U Z T Z 



2 

In this case, the quantity analogous to (|58l), which tends to a constant at large times, is the ratio 

- 2 ~--^^=- (61) 

V\ V\ 2y^(0) 

Its limiting value depends on the initial conditions. The ratio ( |5"8| ) goes to zero as r -1 with 
asymptotics depending on the initial data. 



4. Landau- Lifshitz equation 



Now we investigate the planar motion of a charged particle in a constant homogeneous elec- 



tromagnetic field in the framework of the so-called Landau-Lifshitz equation [20]. This is an ap- 
proximate equation describing the physical solutions to the LD equation. It is obtained from the 
LD equation by the reduction of order procedure, while A being assumed to be a small parameter. 
Thus, let us seek for solutions to the LD equation (|30| ) in a class of functions which 



d k m a 



dr k 



O(l), 



d k V r 



O(l), k = 0,oo, 



(62) 



with respect to the small parameter A. Also, we restrict ourself by the first correction in A to the 
Lorentz equation. 



15 



Differentiating the first equation in (30) with respect to r, we find m a . Then we substitute it 
to the initial equation and come to 



m. 



-u(£ a - pv a - XpV a - Xpv a )- 



(63) 



By the use of this relation, the second equation in ( p0[ ) describing the evolution of the momentum 
v n has the form 



V a = UEap^V 1 + \u?p{i a - pv a ), 



where we neglect the higher orders in A. Convolving this equation with £ Q , we arrive at 



u = -2\uj z u(£ z + u 



(64) 



(65) 



whence 



(£ 2 + u(0))e 2Al ^ v -u(0)' 
'1,0,0), the substitution ( |36|) gives 



r = ±i; 



u 



«(o) 



l + 2Aw 2 n(0)r : 



e = 0. 



(66) 



If = —U). 



(67) 



If £ a = (0,0, 1), the substitution @ results in 



ip = U). 



(68) 



These formulas are in agreement with Eqs. (47) and (49) up to the first order in A. In the case 
f* = (1,0,1), we have 

VRo) 



Vq - V 2 



i>i(0) + y/vffyuT 
a/1 + 2Aw 2 u(0)r ' 
_ v Q (0) - v 2 (0) + 2w(ui(0) + Awv^(0))t + ^l(0)w 2 r 2 



^l + 2Aw 2 u(0)r' 



(69) 



y^l + 2Aw 2 u(0)r 

We see that the system passes to the universal regime at sufficiently large proper times 



.1/2 



.-1/2 



Vl 



V o + V 2 ~ W 



T 3/2 



^0 - ^2 ~ W 



(70) 



'2A' ~ V2A ' " V2A 

The limiting value of the ratio (pq) is the same for these solutions as for the physical solutions to 
the exact LD equation. Thus, the Landau-Lifshitz equation correctly reproduces the asymptotics 
of the physical solution at large r in spite of the fact that this asymptotics is not regular in A. 
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IV. DISCUSSION 



In this paper, we have investigated a planar motion of charged particles obeying the LD equation. 
We gave a detailed study of asymptotics of the planar physical solutions to the LD equation in a 
constant homogeneous electromagnetic field. One of the main results of our study is an interesting 
asymptotics of the physical solution to the LD equation in this electromagnetic field with vanishing 
invariants I\ = I2 = 0. According to the classical radiation reaction theory, the charged particles 
moving on the plane in such the electromagnetic field must have the same ratio ( |58| ) of the momenta 
components at sufficiently large times. The existence of this asymptotics can be verified in a purely 
quantum electrodynamical context. 

Namely, the Dirac equation in that electromagnetic field can be exactly solved pS|| . Then, 
by the use of this complete set of solutions, we construct quantum electrodynamics on the given 
background p9| ] and define the S-matrix. The operators of covariant momenta Vq = mxiQ and 
V\ = mvi commute. Therefore, we can construct a complete set of their eigenfunctions and 
evaluate the transition amplitudes to these states. For the one-photon radiation amplitude, the 
transition probability reads as 

Y,m ] \Vo,V 1] k,\){V Q ,V 1] k,\\U\p), (71) 

k,\ 

where U is the evolution operator over an infinite time, f3 are the quantum numbers characterizing 
the initial state of the electron, and k and A are the momentum and polarization of the radiated 
photon. Provided the classical radiation reaction theory is viable, there must exist a range of the 
quantum numbers j3 and strengths of the electromagnetic field such that the transition probability 



(71) is mostly concentrated on the cubic parabola 

2e 2 E 



V a = --^Vl (72) 



Of course, using formula (71), we disregard the multiple photon production and the production of 



the electron-position pairs, but these contributions are proportional to the fine structure constant 
and neglible for reasonable strengths of the electromagnetic field. We postpone a detail study of 
this effect in quantum electrodynamics for a future research. 



Also note that, as it was pointed out in [13|, the notion of a physical solution seems to be quite 
general and can be exploited to give a proper interpretation of the higher-order derivative corrections 
to effective actions in quantum field theory (see, e.g., [30-3^] and also their quantization in 1 33 , |34|). 



Usually, these terms arise from the heat kernel expansion over the regularization parameter A or 
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large mass m of the one-loop correction |35|| whereas the higher terms of this expansion are casted 
out. By analogy with the considerations presented in Appendix, we should demand a regularity 
of solutions to the effective equations of motion in the small expansion parameter - the inverse 
regularization parameter or inverse large mass. Then a neglect of the nonlocal reminder of the 
effective action is justified. In many instances, this regularity can be related to the regularity with 
respect to the coupling constants. Though, it does not mean that we assume a smallness of the 
couplings or these higher-order derivative terms. As far as the heat kernel expansion is concerned, 
one can distinguish the corrections of three types: i) the divergent higher derivative terms like, for 
example, in R 2 gravity; ii) the higher derivative terms disappearing in the regularization removal 
limit; iii) the higher derivative terms resulting from the large mass expansion of a finite part of the 
effective action. In these cases, the coefficients at the higher derivative terms have a form 

i)~X^ + e a C\)/fa(m/A), ii) e 6 (A)/ 6 (m/A), iii) e c {\)f c {R/m 2 ), (73) 

where e and / are some functions going to zero at the origin, R schematically denotes the field 
strengths or their derivatives. The coefficients in the first case are the renormalized inverse couplings 
A~ . Now it is easy to see that a regularity of the expression in the coupling constants A a implies its 
regularity in A -1 . As for the large mass expansion, we additionally have to require a regularity of 
solutions to the effective equations of motion in m _1 . Just to demonstrate what we mean, consider 
the effective action with a higher derivative correction coming from the (self) interaction 

S[(f)] = ^J dxcj)(-n -m 2 + \\3 2 )0, (74) 

where A is proportional to the coupling constants or, possibly, to the inverse large mass. The 
physical sector of this model is equivalent to 

S phys [<f>] = \j ' dx<t> [-□ - (v / T+4W - 1)/2A] 0. (75) 

As long as as the coupling constants are scalars with respect to a symmetry group of the model, 
its physical sector possesses this symmetry as well. An elimination of the unphysical sector in 
free models is a simple task, but, for full interacting models, such an explicit elimination becomes 
complicated and results in an appearance of the nonlocal projectors to the physical states in the 
effective action. 



Appendix A: Regularity condition in Maxwell electrodynamics 



In this appendix, we show that the solutions to the coupled system of Maxwell-Lorentz equations 
are regular in the coupling constants before taking a regularization removal limit. 
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Consider a model with the action functional (g). The effective equations of motion of a charged 
particle with a bare mass fh look like 



mx,(T) = eF t , u (x(T))x»(T) + 4e 2 J dr'e(X (T,r'))5' \X 2 (r,r') - s)X [ ^t,t')x v] (t')x"(t), (Al) 

where X^(t,t') : = x^(r) — x„(r'), and we have used the regularization of the retarded Green 
function of the form 

G-(x) = ^(z 2 ) -> G~(x) = ^(x 2 -e), (A2) 

where e is a regularization parameter. We shall assume that F /11/ (x) is smooth and bounded on the 
spacetime, F fJil/ (x(T)) vanishes at r < To for some to, and the solution to Eq. (Al) has the following 
asymptotics in the past 

Xfiir) = (t - t )v^ + x^(t), v 2 = l, (A3) 

where x^(t) is zero at t < To, and is a constant 4- vector. The unphysical solutions may appear 
on the scales of a classical electron radius. Therefore, we are looking for the solution of the form 

x fj.( T ) = r e y^,(T/f e ), f e := e 2 /m. (A4) 

The function y^(s) also depends on other dimensionless combinations of parameters entering Eq. 



(Al), but we do not write them explicitly. Upon substitution (|A4|), we arrive at 

Us) = ^f e F liv (y( S ))f(s) + 4 [ ds'9(Y°(s, s'))5'(Y 2 (s, s') — e)Y[ M (s, s')y v }(s')y v (s), (A5) 



where e := e/r 2 and other evident redefinitions have been done. The integral in this equation can 
be taken 



Y { ^s,s')y u] (s, S ') Y bl (s, S ')y v] ( S ') 

(y p (s')YP(s, O) 2 + (y p (s')YP(s, s')) 3 [ M ' [ ' ' { 



r(s), (A6) 



where s' is determined by the conditions Y 2 (s, s') = e, Y°(s, s') < 0. The LD equation is obtained 
by expanding Eqs. ( |A5| ) and ( |A6| ) in the asymptotic series in e around zero and discarding the terms 
vanishing at e — Y 0. In view of the asymptotic behaviour of the solution ( |A3|) , the asymptotics of 
the integral I at small f e readily follows 

I _^e-V 2 fly {ll {s)v v] y v {s). (A7) 

Hence, in this limit, the integro-differential equation ( |A1[ ) reduces the Lorentz-type differential 
equation with the effective strength tensor: F^ v plus the correction ( |A7| ). According to the general 
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theorems of ordinary differential equations theory, the solutions to this equation are regular in f e . 
The effective strength tensor tends to zero when f e — > 0, and, in this limit, the particle moves along 
a straight line 

y^s) = 0. (A8) 

Restituting the notation (|A4|), we see that the limiting trajectory of the particle is indeed regular 
in f e . If we compel the Lorentz force to be constant at f e — > increasing the strength of the 
electromagnetic field, the equation of motion ( |A5| ) turns into the ordinary Lorentz equation and 
its solutions cease to depend on r e . Another possibility is not to scale with f e as in (A4), but 



consider Eq. (Al) with the solutions of the form x^(r /r e ). Drawing the same reasonings, it is easy 
to see that these solutions are regular in f e as well. As long as 

Te ~ 1 + bfe/e 1 / 2 ~ 1 + b- x e x l 2 /f e ' ( } 

where b is some constant, the regularity of a solution in f e implies its regularity in the renormalized 
classical electron radius r e and the regularization parameter e 1 / 2 . 
If we used another regularization of the Green function 

2tte 

then the asymptotics of the integral / would become 

/>oo 

/ _> Ae-^fly [lA (s) v v] y"(s) / dtg'{t 2 ). (All) 

The above arguments applied to this case reveal that there exists a regular limit r e — > of solutions 
to the effective equations of motion of a charged particle with the regularization ( A10| ) as well. 



G~(x) = d -^9(x')g(x'/e), / ,l.n,(.r) I. (AH)) 



A regularity of the electromagnetic field generated by this charged particle is also obvious. The 
regularization of the Green function is equivalent to the regularization of the current 

f{x) -> f{x) = U x j dyGfix - y)fiy), (A12) 

where ja{x) is the current of a point charge. The regularization Gfix) of the advanced Green 
function is analogous to the retarded one ( A10| ). The effective equations of motion of a charge with 



the regularized current have the form (Al), but with the "effective" Green function (see for details, 
e.g., @) 

G e J f ix -y):= J dzdz'U x G- £ ix - z)G~{z - z')U zl G+(z' - y). (A13) 
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It is not difficult to show that the effective regularized Green function is zero in the past light cone. 
A Poincare-invariance of this Green function implies 



G e J f (x) = e -^-e{x 2 )~g{x 2 /e) + e- x e(-x 2 )h(-x 2 le), (A14) 



where g(x) is a generalized function satisfying the condition ( |A10 ), while 



dxh(x) = 0. (A15) 



o 



The part of the effective Green function ( |A14|) with the support lying outside of the light cone 
does not contribute to the integral in Eq. (|A1|). Thus, we revert to the case of the Green function 
regularization considered above. Notice also that we can use the retarded Green function of the 
form (|A^) in the regularized current ( |A12 ) instead of the advanced Green function. It can be proven 



that the effective Green function takes the form ( A14| ) in this case too 



When we pass from Eq. ( |A1| ) to the LD equation (||), we turn to a "truncated" description of 
a charge. Because of the truncation, spurious solutions arise which are non-regular in the coupling 
constants. They should be excluded, since only the regular solutions, which we call physical, to the 



LD equation are close to the solutions of (Al) at small the regularization parameter e. 

Indeed, the LD equation is derived from ( |A1[ ) breaking off the series in e 1 / 2 . A regular in r e 
solution is regular in e 1 / 2 . If we substitute such a solution of the LD equation to Eq. ( |A1| ) expanded 
in the series in e 1 / 2 , we shall ascertain that the obtained expression tends to zero with e — > for any 
r. On the other hand, if we substitute a non-regular solution of the LD equation to Eq. ( |A1| ) then 
the limit e — > of the obtained expression may not exist. Making a general non-regular in r e ansatz 
to the LD equation (JSj), it is not difficult to see that the non-regular solutions to LD equation have 
to possess an essentially singular point at r e = provided some miraculous cancelations do not 
occur. Since it is the LD force which is responsible for the singularity, it should dominate over the 
Lorentz force at small r e and certain proper times r. Therefore, the non-regular solution to the LD 
equation tends to the non-regular solution ( |l4| ) to the free LD equation at r e — > 0. It can be directly 
verified from (||) stretching the proper time r — > r e r. The solution ( |l4| ) does own the essentially 
singular point at r e = 0. But if we substitute this solution to Eq. (|Al|), it blows up at e — > (the 
regularization parameter also enters the renormalized mass). So, the non-regular solutions cannot 



be regarded as approximate solutions to (Al) at small the regularization parameter e. 
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